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Abstract. Anderson localization has been a subject of intense studies for many years. In this context, 
we study numerically the influence of long-range correlated disorder on the localization behavior in one 
dimensional systems. We investigate the localization length and the density of states and compare our 
numerical results with analytical predictions. Specifically, we find two distinct characteristic behaviors in 
the vicinity of the band center and at the unperturbed band edge, respectively. Furthermore we address 
the effect of the intrinsic short-range correlations. 
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1 Introduction 

Over the past fifty years the Anderson model of local- 
ization has become a paradigm for investigations of elec- 
tronic transport in the presence of static disorder [l][2]. 
Depending on the strength of the disorder the wave func- 
tions of non-interacting electrons vary from delocalized 
to exponentially localized. This so called Anderson local- 
ization was demonstrated recently in experiments with 
matter waves of Bose-Einstein condensates in random or 
quasi-periodic optical potentials 'S'.'i'. Due to significant 
electron-electron and electron-phonon interactions in real- 
istic materials, e.g. imperfect crystals, the direct observa- 
tion of Anderson localization is rather difficult. Therefore 
the success of the model is driven to a large extent by the 
universal properties of the resulting phenomena. 

The characteristics of the disorder potential have a 
prominent influence on Anderson localization [2] . For one 
dimensional (ID) systems with a spatially uncorrelated 
potential it has been shown that all electronic states are 
exponentially localized [5]. Later it was realized that by 
introducing correlations, one can partly suppress the lo- 
calization, at least for weak disorder [6]. Moreover, there 
have been discussions of a delocalization-localization tran- 
sition even in ID for long-range correlated disorder poten- 
tials [7||9]. However, it has early been recognized that the 
apparent transition results from a rescaling of the disor- 
der potential by normalizing the variance, which becomes 
system-size dependent for the considered self-affine poten- 
tial landscapes 
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In the present work, we focus on long-range power law 
correlated potentials with a correlation function C{£) cx 
and correlation exponents a > 0. The variance of the 
resulting potential is independent of the system size. This 
so called scale-free disorder is by no means artificial. It 



appears naturally in a large variety of physical systems 

The outline of this article is as follows. In the next 
section we briefly review the main results for Anderson 
localization in ID and in the presence of correlated disor- 
der. We also discuss the numerical methods used in our 
calculations. In Sec.|3]we present and discuss the obtained 
results. Finally, the article is summarized in Sec. |4j 



2 Model and Numerical Techniques 

2.1 Anderson Model of Localization with Long-Range 
Correlated Disorder 



The Anderson model [I][2] is widely used to investigate the 
phenomenon of localization in disordered materials. It is 
based upon a tight-binding Hamiltonian in site represen- 
tation, 

where \i) is a localized state at lattice site i and {ij) de- 
notes a restriction of the sum to i and j being nearest 
neighbors. The hopping parameter t = It? / {2m* a?) is de- 
fined in terms of the lattice spacing a and the effective 
mass m* . In the following, we set t = 1 and thus fix the 
unit of energy; lengths are measured in units of a. The 
on-site potentials £i are random numbers, chosen accord- 
ing to some probability distribution P(e) characterized by 
the mean (e^) and the correlation function {eiEi+i) = C{£). 
However, usually the site energies are taken to be statisti- 
cally independent. For example, a common choice of -P(e) 
is a Gaussian white noise distribution or a box distribu- 
tion of width W, both with (e,) = and C{e) = ^Se^. 
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Other distributions have also been considered [2 14 ■ 16 
In the present work we are interested in the innuence of 
long-range correlated disorder potentials on the Anderson 
model, i.e., using a correlation function of the form 



(2) 



where a is the correlation exponent which determines the 
strength of the correlations. The associated spectral den- 
sity is given by 



OC \q\ 



(3) 



In order to generate physically reasonable long-ranged cor- 
related potentials, the correlation function ([2]) should de- 
cay with distance I and therefore a > 0. For ID sys- 
tems the correlations are considered to be long ranged 
for < a < 1. The value a = 1 corresponds to the case 
of uncorrelated disorder. Notice that the delocalization- 
localization transition observed in Refs. [7]-[9] relies on 
power-law correlations with a < 0. These correlations in- 
crease with distance i and lead to a system-size depen- 
dence of the disorder potential. Thus a rescaling of the 
results is required which yields the corrected localization 
behavior without transition 
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In general, it is extremely complicated to obtain an- 
alytical results of transport properties for the Anderson 
model of localization. For example, only in the ID case 
rigorous proofs of strong localization for all energies and 
disorder strengths have been given [5|. Moreover, the ex- 
plicit energy and disorder dependence of the localization 
length A for weak disorder has been derived 17p8 . There 
are also some results for ID systems with long-range cor- 
related disorder. For energies close to the band center 



^ 2) a weak disorder expansion for A yields 



A{q,W) 



\sin\q) 12 



S{2q) 



(4) 



where E{q) = — 2cos(g). It is important to notice that 
even for the correlated case the localization length remains 
proportional to W~^. Thus, for small disorder strength 
and in the vicinity of the band center, the presence of 
correlations leads to an effective disorder strength WcS = 
W ^J~S{2q)^ which depends on energy via the spectral den- 
sity S. 

On the other hand, close to the unperturbed band edge 
2) a different universal behavior can be found 



\E\ 



19 



20 . Here, the density of states (DOS) and the localization 
length are characterized by units of energy e and length 
A, 



A^{E) = Xgc.{E/e) , 



(5a) 
(5b) 



where /„ and ga are universal functions. For uncorrelated 
potentials, i.e. a = 1, there is a closed analytic form avail- 
able for both functions 21 . The units e and A can be 



expressed in terms of the hopping parameter t and the 
disorder strength w = W/y/l2, 



A 



y4/(4-a) 
„-2/(4-a) 



(6a) 
(6b) 



In particular, it follows that at the unperturbed band edge 
(1^^1=2), 



AU\E\ 
Pc.i\E\ 



2,w) OC w-y , 
2,W) ocw-y 



(7a) 
(7b) 



with y = 2/3 for a = 1 and y = 2/(4-a) for a < 1 [2Tp2] . 
Most interestingly, within the so called white-noise moHel 
(WNM) scaling expressions identical to Eqs. ^ are ob- 
tained where a is replaced by the dimension d of the sys- 
tem. We would like to emphasize that the WNM does not 
consider long-range correlated disorder, but the similarity 
suggests that the influence of power-law correlations at 
the band edge may possibly be interpreted in terms of an 
effectively reduced dimensionality. 



2.2 Numerical Methods 

In order to generate random numbers with power-law cor- 
relations ([2]) we use the modified Fourier filtering method 
(FFM) 23 with one additional step. Hereby, the correla- 
tion function is given by 



C{£) = (l+f) OC 



[i » 1) 



(8) 



which avoids the singularity at = and therefore im- 
proves the filtering. For large distances £ ^ 1 this form 
resembles the desired power-law behavior ([2]). An addi- 
tional benefit of Eq. ([s]) consists in avoiding aliasing ef- 
fects, which may obscure the long-range character of the 
correlations. The spectral density Sicf) can be calculated 
analytically and is given in terms of a modified Bessel 
function [23]. Overall, the special choice of the correlation 
function ffST leads to a modified behavior with respect to 
the correlation exponent a. While for the pure power-law 
the case of uncorrelated disorder is obtained for a = 1 (see 
also Eq. ([3])), the modified FFM generates uncorrelated 
random numbers only in the limit a — ^ 00. Therefore, it is 
expected to find an infiuence on the localization behavior 
even for a > 1. 

Additionally, we shift and normalize the obtained se- 
quence of correlated random numbers such that the mean 
vanishes and the variance is W'^/12 ^9j. Thereby, only the 
strength of the disorder is adjusted while the correlations 
of the potential are preserved. Then we calculate the lo- 
calization length using a standard transfer-matrix method 
(TMM) U with a new seed of the random number gener- 
ator for each parameter combination {E, W, a). The DOS 
is calculated by diagonalizing the Hamiltonian ([l]) and 
counting the eigenvalues in finite energy bins. In each case 
100 realizations are taken into account for averaging. 

Due to the symmetry of the Hamiltionian ([T]) and the 
chosen symmetric disorder distribution P{e) the results 
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Fig. 1. (Color online) Localization length Aa vs energy E at 
W = 0.5 for a ID chain of length L = 2^^ obtained from TMM 
calculations (a). Full lines correspond to Aa given by Eq. Q. 
DOS pc for a ID chain of length L = 2" (b). The dashed line 
denotes the exact DOS for a ID system without any disorder, 
p — 1 1 (7rV4 — -E^). Different symbols/colors denote various 
correlation parameters a. 

should depend on the absolute value of E only. This is 
confirmed by our calculations. Therefore we use the arith- 
metic average of values for —E and ^E to improve the 
accuracy of our results. 

In the following, we consider chains of fixed length 
L = for the TMM calculations and L = 2^^ for the 
DOS computation. The correlation exponents are a = 
00,2.0,1.0 and 0.5. 



3 Results and Discussion 

In Fig. [T] we show the energy dependence of A and the 
DOS for fixed disorder strength W = 0.5 and different cor- 
relation strengths a. For energies inside the unperturbed 
band {\E\ < 2) the localization length is increased by long- 
range correlations. In this region one observes an overall 
good agreement with the weak-disorder result given by 
Eq. Q. On the other hand, near the band edge {\E\ ~ 2) 
A is decreased by correlations, which will be discussed in 
more detail below. Another important quantity for trans- 
port studies in general is the DOS p{E). For an ordered 
system the DOS can be calculated analytically. In this 
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Fig. 2. (Color online) Difference of the DOS pa for a ID chain 
of length L = 2^^ and W — 0.5 compared to the uncorre- 
lated case. Different symbols/colors denote various correlation 
parameters a as given in the legend. 



case it develops van Hove singularities at the band edges 
\E\ = 2, which is typical for ID systems [24] and indi- 
cates the presence of long-range order [2]. In case of un- 
correlated on-site energies the van Hove singularities are 
smeared out while maintaining peaks at \E\ = 2. Inside 
the unperturbed band the DOS is only weakly changed 
compared to the DOS of the ordered system. For \E\ > 2 
there is an exponential band tail [2]. For correlated on- 
site energies Fig. [2] illustrates that the van Hove peak is 
even less pronounced and shifted towards the band center, 
which results in a decreased DOS at \E\ = 2 with decreas- 
ing a. On the other hand, in the band tails the DOS is 
getting larger for decreasing a. 

In the following we discuss the two regions separately 
and in more detail. Firstly, we focus specifically on the 
behavior in the center of the band. Figure [3] shows the 
disorder strength dependence of the localization length at 
E = 0. One can see that introducing correlations leads to 
a systematically larger localization length. Only for very 
strong disorder the influence of the correlations vanishes 
and the localization length for finite a approaches the un- 
correlated localization length. Also shown is the localiza- 
tion length obtained from Eq. Q using the same parame- 
ters as in the TMM calculation. It is important to notice, 
that there is no fitting procedure involved at this point. 
Therefore, Eq. Q allows for an independent check of the 
numerics and in particular of the correct behavior of the 
generated correlations. Figure |3] shows a good agreement 
of A and the weak disorder expansion result Q for 11/ < 1. 
In contrast to A the DOS in the band center has no par- 
ticular dependence on the correlation strength for weak 
disorder as demonstrated in Fig. [2] Therefore the infiu- 
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Fig. 3. (Color online) Localization length Aa vs disorder 
strength W at the band center E = Q for a ID chain of length 
L = 2^^ obtained from TMM calculations. Different symbols 
denote various correlation parameters a. Dashed lines show 
given by Eq. Q. 



ence of the long-ranged correlations can be interpreted in 
accordance with Eq. (|4| as an effective disorder of the 
potential Hereby, the effective disorder strength becomes 
smaller with decreasing a. 

A quahtatively different behavior of the localization 
length is observed for energies close to the unperturbed 
band edge (|£'| — 2). This is shown in Fig. l4] Here the 
locahzation is enhanced in case of a correlated potential 
and A is consequently smaller compared to the uncorre- 
lated case. However, for weak disorder W <1 one can still 
observe a power-law decay according to Eq. Q , but with 
an exponent y < 2/3. The results of least-squares fits of 
Eq. Q to the numerical data are summarized in Tab. [l] 
and also shown in Fig. |4j The obtained exponents for the 
uncorrelated potential as well as for a — 2.0 and 0.5 agree 
with the values predicted by Eqs. ([T]). Only for a = 1.0 
we observe a distinct deviation from the expected value 
y = 2/3 which could be a consequence of the modified 
FFM [23]. From Eq. ^ one expects the FFM to generate 
an uncorrelated potential for a = 1 since S{q) becomes 
independent of q in this case. However, due to the short- 
range cutoff introduced in Eq. ([s]), the spectral density is 
given in terms of a modified Bessel function [231, which 



is independent of q only in the limit a — >■ oo. Therefore, 
the behavior of the localization length and the DOS at the 
band edge is sensitive to the specific choice of the short- 
range cutoff for sufficiently weak disorder. On the other 
hand, for very strong disorder the localization lengths for 
the correlated cases coincide with the uncorrelated case. 
Figure [5] shows the DOS at the band edge for small dis- 
order strengths. Similar to the localization length one ob- 
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Fig. 4. (Color online) Localization length Aa vs disorder 
strength W at the unperturbed band edge E = 2 for a ID 
chain of length L = 2^®. Different symbols indicate results ob- 
tained from TMM calculations. Error bars are smaller than the 
symbol size. The full line shows the exact result for uncorre- 
lated disorder [2l]. Different colors denote various correlation 
parameters a. Dashed lines show results of a least squares fit 
of Aa oc W~''' to the data. 



serves a power-law decay of the DOS, which is also in 
accordance with the scaling law The fitted exponents 
are listed in Tab.[T] Their values agree with the results ob- 
tained from the localization length but are less accurate. 



4 Summary and Conclusions 

In summary, we have numerically investigated the role of 
long-range correlated disorder on Anderson localization in 
ID systems. Specifically, we have studied the behavior of 
the localization length and the DOS as functions of en- 
ergy and disorder strength for different correlation expo- 
nents. Hereby, we found two regions in the electronic band, 
where universal behavior can be observed. In the vicinity 
of the band center and for weak disorder the localization 



Table 1. Results of a least squares fit of Aa oc W~^^ and 
Pa oc W~^'' to the numerical data for weak disorder. For Aa 
and Pa the range W G [0.1,0.5] and W € [0.4,1.0] was used, 
respectively. The obtained exponents are compared to the ex- 
pected values y according to Eqs. ([7|. 



a 


oo 


2 


1 


0.5 


y 


2/3 « 0.67 


2/3 ^ 0.67 


2/3 « 0.67 


4/7 « 0.57 


VA 


0.68 ±0.01 


0.66 ±0.01 


0.60 ±0.02 


0.55 ±0.01 


Vp 


0.65 ±0.01 


0.61 ±0.03 


0.56 ±0.02 


0.54 ±0.03 
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Fig. 5. (Color online) DOS pa vs disorder strength W at the 
unperturbed band edge E — 2 for a ID chain of length L — 2^^. 
Different symbols indicate various correlation parameters a. 
The full line shows the exact result for uncorrelated disorder 
|21| . Dashed lines denote results of a least squares fit of pa oc 
to the data. 



length is determined by an effective disorder strength. The 
localization length remains proportional to W~'^ in this 
regime. The explicit expression, which has been derived 
in Ref. [6], is given by Eq. Q. A qualitatively different 
behavior is found at the band edge. Here, the localiza- 
tion length and DOS are given by Eqs. ([5]), which implies 
a power-law behavior for both quantities as functions of 
disorder strength. A comparison with the WNM suggests 
that in this region the long-range correlations lead to an 
effectively reduced dimension. 

Overall, we find good agreement of our numerical re- 
sults with the respective analytical expressions. However, 
we also find that the particular realization of the long- 
range correlations plays an important role. In the present 
case the FFM 
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makes use of a modified correlation 
function given by Eq. ([8|, which contains the desired long- 
range but also short-range contributions. Our numerical 
results indicate that the modification leads to observing 
an unexpected infiuence of the correlations even for cor- 
relation exponents a > 1. This behavior should be taken 
into account for investigations involving realistic disorder 
potentials with long-range correlations. 
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